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Abstract 

The fragmentation functions D a ^ff(x, /J* 2 ) of a heavy hadron H, with mass run, 
satisfy the phase-space constraint D a ^H(x, /x 2 ) = for x < m 2 H /ii 2 , which is vi- 
olated by the naive /i 2 evolution equations. Using appropriately generalized /j, 2 
evolution equations, we reconsider the inclusive hadroproduction of prompt J ftp 
mesons with high transverse momenta in the framework of the factorization formal- 
ism of nonrelativistic quantum chromo dynamics, and determine the resulting shifts 
in the values of the leading colour-octet matrix elements, which are fitted to data 
from the Fermilab Tevatron. 

PACS numbers: 13.85.-t, 13.85.Ni, 13.87.Fh, 14.40.Gx 



*Permanent address: II. Institut fur Theoretische Physik, Universitat Hamburg, Luruper Chaussee 
149, 22761 Hamburg, Germany. 



In the framework of the Q CD-improved parton model, the inclusive production of 
single hadrons is described by using fragmentation functions, D a _+ h (x, /i 2 ). The value 
of D a ^h{x, p 2 ) corresponds to the probability for a parton a which comes out of the 
hard-scattering process to form a jet which contains a hadron h carrying the longitudinal- 
momentum fraction x = (Ph+Ph) I (Pa+PDi where p a and ph are the four-momenta of a and 
h in the infinite-momentum frame. Here, p is the fragmentation scale, which is typically 
chosen to be of the order of the centre-of-mass energy sfs (transverse momentum of h) 
in lepton-lepton (lepton-hadron and hadron-hadron) collisions. At next-to-leading order 
(NLO), p is identified with the mass scale at which the collinear singularities associated 
with the outgoing parton a are factorized. In the case of a light hadron, with mass 
rrih <C p, the p 2 dependence of D a ^h(x, p 2 ) is determined by the timelike Altarelli-Parisi 
evolution equations, 

^D a ^ h (x,p 2 ) = E f -P^ b {y^ 2 )D b ^ h ( X -,A , (l) 

where P^ b (x, p ) are the timelike a — > b splitting functions. Ready-to-use expressions 
for P^_; h (x, p 2 ) through NLO are collected in the Appendix of Ref. HI] 

In Eq. ([I]), x may, in principle, be arbitrarily low for any value of p. In the case of a 
heavy hadron H, with mass rrin^p, this conflicts with the fact that the fragmentation 
process is only allowed by kinematics if the virtuality p 2 of the fragmenting parton a 
satisfies the condition p 2 > m 2 H /x. This may be understood as follows. Consider the 
general situation where a fragments into n particles, with masses and four-momenta 
Pi, and choose the coordinate system so that p^ = 0, 0,p^). Defining the longitudinal- 
momentum fractions as Xj = +pf)/(Pa +pI), we then have 

vl = ( P ° a + pI) E (rf - Pl) = E (2) 

i=l i=l X i 



where pT,i = y (pj) + (pf) are the intrinsic transverse momenta. If all final-state par- 
ticles, except for H, are light and the intrinsic transverse momenta are neglected, as is 
usually done in the parton model, then it follows that p\ > m 2 H /x. If we identify p 2 = p 2 , 
then this phase-space constraint leads to the condition D a ^ H (x, p 2 ) = for x < m 2 H / p 2 . 

In order to properly implement this condition, Eq. ([l|) must be generalized. Such a 
generalization was proposed in Refs. 0, |3| for the p 2 evolution of D g ^ H (x, p 2 ) generated 
by the g — > g splitting. It is straightforward to extend this formalism to also include the 
fragmentation of the other partons and the nondiagonal evolution effects. This leads us 
to the ansatz 

D a ^ H (x,p 2 ) =E I n ~t [ —G a ^ b (y,q 2 ;p 2 )d b ^ H f-,<? 2 ) , (3) 

1 There is an obvious typographical error in the published version of Ref. [Q, which was absent in the 
preprint version thereof. In the line before the last of Eq. (17), lnln(l — x) should be replaced by ln(l — x). 
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where G a ^b(y, q 2 ', ^ 2 , ) corresponds to the probability for the parton a, with virtuality 
/i 2 , to emit a parton b, with longitudinal- momentum fraction y and virtuality q 2 , and 
db^H^z, q 2 ) to the one for parton b to subsequently decay to a hadron H, with longitudinal- 
momentum fraction z relative to b. The fi 2 dependence of G a ^b(y, q 2 ] fJ> 2 ) is determined 
by the evolution equation 



^ 2 G a -,b(y, g 2 ; /i 2 ) = E Jl ^ T Jc(z, f)G c -, b (| , g 2 ; V) , (4) 



which is similar to Eq. ([[J), except that the virtuality of the intermediate parton c is taken 
to be zfi 2 instead of /i 2 .[] Furthermore, G a ^b(y, q 2 ] l^ 2 ) satisfies the boundary condition 

G a -+ b (y, /j, 2 ; /i 2 ) = 5 ab 5(l - y). (5) 

According to the above argument, G a ~, b (y, q 2 \ ^ 2 ) is subject to the phase-space constraint 
G a ^b{y,q 2 ] fi 2 ) = for y < q 2 /fi 2 . For the same reason, we have db^H(z,q 2 ) = for 
z < m 2 H /q 2 . From Eq. (||) it hence follows that D a ^n(x, fi 2 ) = for x = yz < m 2 H /fi 2 as 
desired. 

In the perturbative calculation of the fragmentation function D a ^n (x, /1q) at the initial 
scale fj,Q, with fio^mn, <i a ^#(x, g 2 ) acts as a source density, in the sense that 



_,2 

D a ^ H (x,n 2 ) = I " ^-d a ^ H (x,q 2 ), (6) 



.2\ r° dq 2 , , 



which follows from Eq. @ by approximating G a ^b(y , Q 2 ]^ 2 ) with the aid of Eq. (|5|). In 
practice, the upper bound of integration in Eq. (||) is taken to be infinity because the 
integrand rapidly falls off with increasing q 2 M, [51] . 

Differentiating Eq. (Q) with respect to ln/x 2 and substituting Eqs. (£|) and fl5|) on 
the right-hand side, we may eliminate G a ^b(y, Q 2 ', ^ 2 ) and thus obtain a single set of 
inhomogenious integro-differential evolution equations for _D a _^#(x, /i 2 ), namely, 

d ;D a ^ H (x,ii 2 ) = d a ^ H (x,/j, 2 ) +J2 I — p S-lb{y^ 2 ) D b-^H (-,2/Ai 2 ) , (7) 



which is to be solved imposing the boundary condition D a ^n {x, Tnjj) = 0. Evidently, the 
solutions of Eq. (|7|) satisfy the condition D a ^n(x, jj 2 ) = for x < m 2 H / n 2 as they should. 
Notice that Eq. (0) carries over to NLO as it stands. One just needs to include the NLO 
corrections to d a ^ H (x, /i 2 ) and P^} b (y , /i 2 ) and to employ the two-loop formula for the 
strong coupling constant a s ([i 2 ). 

We now explore the phenomenological implications of the modified fi 2 evolution for 
the case of prompt J/^-meson production via fragmentation, H = J/ip, in the framework 
of the factorization formalism of nonrelativistic QCD (NRQCD) ||. At the starting 
scale fiQ, the leading contributions arise from the fragmentation processes a — > cc[n] 

2 Here, we deviate from Eq. (4) of Ref. M, which refers to the evolution in q 2 . 
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Table 1: Values of the leading J/i[) colour-octet matrix elements (in units of 10~ 3 GeV 3 ) 
and of r resulting from the LO and HO-improved fits to the CDF data [113 with naive 
and modified /x 2 evolution. 





naive 


modified 


LO 


HO 


LO 


HO 




4.58 ±0.73 


2.58 ±0.42 


5.01 ±0.79 


2.86 ±0.46 




104 ± 10 


19.1 ± 2.8 


99.9 ±9.8 


17.0 ±2.7 


r 


3.44 


3.53 


3.43 


3.52 


Xdf 


0.52 


0.52 


0.52 


0.41 



with a = g,c,c and [n] = [1, 3 5i] , [8, 3 Si], where 1 and 8 label colour-singlet and colour- 
octet states, respectively, and the spectroscopic notation 2S+1 L j indicates the spin S, the 
orbital angular momentum L, and the total angular momentum J. The corresponding 
fragmentation functions D a ^j/^ (x, /i 2 ,) may be found in Refs. ||, |jj and their source 
densities d a ^j/^(x,q 2 ) may be gleaned from Ref. M. The fragmentation functions for 
a — u,u, d, d, s, s are generated via the \j? evolution and coincide. Furthermore, we have 
D c ^j/ip(x, /i 2 ) = D 5 _>j/tp(x, /i 2 ), so that we only need to distinguish the three cases a = 
g,c,u. In Figs. |l]a-c, we study the x dependences of D a ^j/^(x, /i 2 ) at jj 2 = 300 GeV 2 for 
a = g,c, u, respectively, comparing the the naive (dashed lines) and modified (solid lines) 
fi 2 evolutions to leading order (LO). We adopt the LO color-singlet matrix element from 
Ref. || and the LO colour-octet ones from Table [l], selecting those which refer to the naive 
/i 2 evolution. The /i 2 evolutions are performed iteratively in x space. The technicalities 
are explained for the naive and modified /i 2 evolutions in Refs. |7], ||, respectively. As is 
well known [|nj, the naive fi 2 evolution equations ([IJ) break down in the low-x region due 
to the presence of large logarithms of 1/x on their right-hand sides. The most dramatic 
consequence of this deficiency is an unphysical divergence in the gluon multiplicity at 
low x. This directly triggers the singular low-x behaviour of D g ^j/^(x, /i 2 ), which is 
exhibited by the dashed line in Fig. |l]a. Via nondiagonal evolution effects, this feature 
also feeds into the quark fragmentation functions, as may be seen from Figs. |I]b and 
c. By contrast, the results based on the modified /x 2 evolution are devoid of such low- 
x divergences. As per construction, they vanish for x < rn 2 j^/ /i 2 = 0.032. On the 
other hand, they merge with the naive evaluation as x approaches unity because the 
impact of the phase-space constraint then fades out. Both features serve as a welcome 
check for our numerical analysis. Further evidence for the physical soundness of the 
modified /i 2 evolution comes from a recent phenomenological analysis of Z-boson decay 



into prompt J /if) mesons via fragmentation [11], where the large logarithms of the fixed- 



order calculation were resummed by using a Monte Carlo cascade model. Except very 
close to the threshold region, the results thus generated were found to nicely agree with 
those obtained from the modified fi 2 evolution. 
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At present, one of the most important applications of the J/^-meson fragmentation 
functions is to describe the inclusive hadropro duct ion of prompt J/ip mesons with high 
transverse momenta in pp collisions at the Fermilab Tevatron, with y/s = 1.8 TeV. In 
fact, fits JS], [12] to the latest data on pp — > J/tp + X taken by the CDF Collaboration 
p| allow one to place stringent constraints on the leading colour-octet matrix elements, 
(0 J /^[8, 3 Si]), (O j ^[8,%]), and (O j /^[8, 3 Pj]), with J = 0,1,2, in the framework of 
the NRQCD factorization formalism ||. Specifically, the data in the upper px range 
are especially sensitive to (0 J ^[8, 3 Si\^, while the \ow-px data essentially fix the linear 
combination 

M m = / J ^[8, %}) + (<D J/ ^[8, 3 P ] 



m 2 



where r is to be chosen in such a way that the superposition of these two channels is 
insensitive to precisely how they are weighted relative to each other. Due to heavy- 
quark spin symmetry, the multiplicity relation (o J ^[8, 3 Pj]j = (2 J + 1) (o J ^[8, 3 P 
is approximately satisfied. Adhering to the analysis of Ref. [§], we now investigate by 
how much the fit values for (p J ^[8, ^x]) and M r are shifted if we pass from the naive 
to the modified jj 2 evolution The philosophy advocated in Ref. M is adopt the fusion 
picture, where the cc bound state is formed within the primary hard-scattering process, 
in the low-p-r regime and the fragmentation picture, where the cc bound state is created 
from a single parton which is close to its mass shell, in the high-pr regime. We first redo 
the LO and higher-order-improved (HO) analyses of Ref. ||, which are based on slightly 
obsolete parton density functions, using the latest LO and NLO MRST sets [13]. We then 
repeat these calculations with the modified /i 2 evolution. The resulting LO and HO values 
for (C W [8, 3 5i]) and are summarized in Table [I] together with the corresponding 

values of x 2 per degree of freedom, Xdf- The LO and HO values for (O j ^[1, 3 Si\^ 
remain the same as in Ref. [§j . We observe that the effect of the modified /i 2 evolution is 
to increase (C J /^[8, 3 Si]) and to decrease M J r ^ . This may be understood from Figs. |l]a- 
c by observing that the use of the modified fi 2 evolution leads to a reduction of the 
fragmentation functions at low and intermediate values of x. This must be compensated 
by an increase of (O j ^[8, 3 S\]j in order to match the high-pr data, which are described 
in the fragmentation picture. On the other hand, the use of this increased value of 
(O j ^[8, 3 Si] \ in the fusion picture must in turn be compensated by a reduction of 
in order to fit the low-pr data. Within each order of perturbation theory, the shifts in 
/(9 J /^[8, 3 Si]) and due to the change of the evolution mode are relatively modest, 

comparable to the errors on these quantities. This may be understood by considering the 
average x values which are probed by pp — > J/ip + X via fragmentation. At the data point 
of highest pt, with px ~ 18 GeV, we have x = 0.73 ± 0.14 in the colour-singlet channel, 
x = 0.90 ± 0.12 in the colour-octet channel, and x = 0.90 ± 0.12 for their combination. 
As is evident from Figs. [l]a-c, the massive-evolution effect is not yet pronounced at such 
high x values. 

In conclusion, the implementation of the phase-space constraint for heavy hadrons in 
the fi 2 evolution leads to a significant reduction of their fragmentation functions at low 
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and intermediate values of the longitudinal-momentum fraction x. However, current de- 
terminations of the leading colour-octet matrix elements for the J / ip meson from Tevatron 
data of pp — > J/ip + X are only modestly affected by this because the x values which are 
typically probed in this reaction are rather high. 
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Figure 1: D a ^j/^(x, /i 2 ) to LO at /i 2 = 300 GeV 2 as a function of x for (a) a = g, (b) 
a = c, and (c) a = u. The results obtained with the modified fj, 2 evolution (solid lines) 
are compared with those obtained with the naive ji 2 evolution (dashed lines). 
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